Giant Casimir effect in fluids in non-equilibrium steady states 
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In this letter we consider the fluctuation induced force exerted between two plates separated by a 
distance L in a fluid with a temperature gradient. We predict that, for a range of distances L, this 
non-equilibrium force is anomalously large compared to other Casimir forces. The physical reason 
is that correlations in a non-equilibrium fluid are generally of longer range than other correlations, 
even than those near an equilibrium critical point. This giant Casimir force is related to a diver- 
gent Burnett coefficient that characterizes an Onsager cross effect between the pressure and the 
temperature gradient. The predicted Casimir force should be detectable with currently available 
experimental techniques. 
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Fluctuation induced forces are common in nature [if. 
The well known prototype of such a force is the Casimir 
force between conducting plates due to quantum fluctu- 
ations of the electromagnetic (EM) field [4] . In this case 
the overall energy scale is set by Planck's constant, h, 
and the force per unit area, or pressure, is 
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where L is the distance between the plates, c the speed 
of light, and where the minus sign indicates an attractive 
force. Other commonly discussed induced forces involve 
thermal fluctuations where the energy scale is set hy k^T, 
where is Boltzmann's constant and T the tempera- 
ture pj]. These thermal forces are important when the 
fluctuations are large and long range. The first instance 
of this type was noticed by Fisher and de Gennes, who 
considered finite-size corrections to the free energy in a 
fluid near a critical point 3i]. One then finds a scale- 
dependent force per unit area, or critical Casimir pres- 
sure, Pc, that is given by Q 
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where Q (x) is a finite-size scaling function with ^ the 
correlation length. One defines a Casimir amplitude, 
A = lim O (x), which for the Ising-like universality class 
may vary from —0.01 to +2 depending on the boundary 
conditions Note that at larger L, \pc\ > |pem|- That 
is, the fluctuations that cause Pc are effectively of longer 
range than those that cause pem- Similar Casimir forces 
have been predicted to exist generally in equilibrium sys- 
tems, when long-range correlations are present due to 
the existence of Goldstein modes. Systems investigated 
include superfluid helium and liquid crystals JJ. 

It has by now been well established that thermal fluc- 
tuations in fluids in non-equilibrium steady states are 



anomalously large and very long range. The most stud- 
ied case is a quiescent fluid in the presence of a uniform 
temperature gradient, VT. Then the non-equilibrium 
contribution to the temperature fluctuations as a func- 
tion of the wave number, k, is given by [S| 
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Here the temperature gradient is taken to be in the z- 
direction, the plates are located at z = 0, L, and fc|| is 
the magnitude of the component of the wave vector k 
in the direction parallel to the plates. In this equation, 
p is the mass density, Dt the thermal diffusivity, and v 
the kinematic viscosity. This result for the temperature 
fluctuations at hydrodynamic length scales was predicted 
a long time ago 0, Isl and has been verified accurately 
by light scattering [l3] and by shadow-graph exper- 
iments ll|. From Eq. (3) one sees that the intensity 
of the temperature fluctuations diverges as fc~^ when 
/c — )■ 0. The purpose of this Letter is to show that, as 
a consequence, these non-equilibrium temperature fluc- 
tuations will cause a Casimir effect that is even more 
significant than the one induced by critical fiuctuations, 
whose intensity only varies as [l^l- Specifically, for 
the scale-dependent non-equilibrium fiuctuation contri- 
bution, pnEj to the pressure we have found 
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Here Cp is the isobaric specific heat capacity, 7 the ra- 
tio of the isobaric and isochoric heat capacities, and a 
the thermal expansion coefficient. Note that for a fixed 
temperature gradient, this Casimir force actually grows 
with increasing L. This anomalous behavior is a refiec- 
tion of the very long spatial correlations in a fiuid in a 
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non-equilibrium state. Experimentally, it may be eas- 
ier to measure the L dependence of the non-equilibrium 
Casimir pressure by fixing the temperature difi^erence, 
AT, between the two plates and varying the distance 
L. If a uniform temperature gradient is present, then 
Vr = 1ST jh^ and for fixed AT the non-equilibrium 
Casimir effect decreases as for larger L. Also note 
that pne can be positive or negative. 

To understand the physical origin of Eq. (4), we start 
with the general observation that a temperature gradient 
can cause normal stresses or pressures, if nonlinear effects 
are taken into account. Specifically, we consider a non- 
linear Onsager cross effect characterized by a non-linear 
kinetic coefficient that we define as knl- On these general 
grounds one expects a non-equilibrium contribution to 
the pressure that is given by 



terms of a non-equilibrium time-correlation function: 
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In standard transport theory, knl is referred to as a 
non-linear Burnett coefficient 13]. It is well known that 
these transport coefficients, which go beyond ordinary 
Navier-Stokes transport coefficients, do not exist due to 
the presence of long-time-tail (LTT) effects [1, 14|. In- 



deed, at this order one expects that knl has a LTT con- 
tribution that diverges linearly with the system size 15, 
[l6j . To account for such a LTT contribution, one should 
consider knl as the sum of a bare contribution, ^^^l' 
sociated with short-range correlations, and a divergent 
contribution ~ L due to long-range correlations with a 
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Substituting knl — '«nl i^**^ -^^l- (^) 



coefficient 

yields a pressure effect associated with short-range cor- 
relations. Since the ratio i^'^i^/ k^-^L will be of the order 
of (j/ L, where ct is a typical intermolecular distance, this 
contribution is small and can be neglected. Also such 
an effect is one of several at molecular scales, includ- 
ing accommodation of the velocity and kinetic energy of 
the molecules with the wall, which do not satisfy uni- 
versal laws and should not be characterized as Casimir 
effects However, substitution of knl — '^nl-^ into 
Eq. (5) yields a genuine non-equilibrium Casimir pres- 
sure due to long-range correlations: 
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Here Ji is a microscopic longitudinal stress current, 
a microscopic energy current, and the angular brack- 
ets denote an average over an A^-particle local equilib- 
rium (LE) distribution function. This time correlation 
function can then be evaluated in what amounts to a 
mode-coupling a ppr oximation by using techniques devel- 
oped previously 3 S]- The chief difference with equi- 



librium mode-coupling theory is that in this case one 
must use hydrodynamic modes appropriate for the non- 
equilibrium steady state considered here. It also means 
that the modes are constructed for a system finite in the 
z-direction. To obtain an explicit analytic result we as- 
sume stress-free boundary conditions for the fluctuations. 
The net result is Eq. (4). One sees explicitly that this 
equation is proportional to the mode-coupling amplitude 
for the heat conductivity times the one for the longitudi- 
nal, i.e. bulk, viscosity [20|. That is, the non-equilibrium 
pressure is an Onsager cross effect with a divergent ^ L 
Burnett kinetic coefficient. 

The second method to derive Eq. (4) starts from a non- 
linear fluctuating pressure, relating it to the temperature 
fluctuations, and then evaluating the temperature fluc- 
tuations by, for example, the method of non-equilibrium 
fluctuating hydrodynamics [2l|. It is convenient to con- 
sider the pressure as a function of a fluctuating energy 
density e + 5e and a fluctuating number density n + 5n 
and write the pressure in terms of the mean values, e, n, 
and their fluctuations (5e, 5n: 
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Taylor expanding this pressure and averaging the result 
leads to a fluctuation renormalization of p{e,n). The 
fluctuations that lead to Eq. (3) arise from a combina- 
tion of entropy and transverse velocity fluctuations with 
vanishing linear pressure fluctuations 0, [2l| . This im- 
plies that to second order in the fluctuations the non- 
equilibrium pressure renormalization is given by 
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Hence, Eq. (4) represents a calculation of the LTT, or 
divergent, part of knl- It thus follows that a measure- 
ment of the non-equilibrium Casimir pressure, pne, will 
not only verify the existence of a novel type of Casimir 
effect, but would also be a direct measure of a divergent 
Burnett coefficient. 

Our main result, Eq. (4), can be derived by two distinct 
approaches. The first approach uses a general statistical- 
mechanical method [l7| to express the normal stress in 
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In deriving Eq. (9) we have used that vanishing lin- 
ear pressure fluctuations implies 5e = —w5n and 6n = 
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—naST. To evaluate the temperature fluctuations in- 
duced by the presence of a temperature gradient one 
should notice that the variation of the temperature 
fluctuations as function of the wave number in Eq. (3) 
causes large finite-size effects when the fluid layer is 
bounded by plates separated by a finite distance L. These 
finite-size effects have been evaluated both for stress-free 



TABLE I: Estimated Casimir pressures 



boundary conditions [22] and for no-slip boundary condi- 
tions [23|. If we assume stress- free boundary conditions, 
for which an explicit analytic expression is available 22] , 
we obtain 
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Upon substituting Eq. (11) into Eq. (9) and using some 
thermodynamic manipulations we find Eq. (4) again. 

To estimate the non-equilibrium Casimir pressure, 
Pne, implied by Eq. (4), we consider as an example liq- 
uid toluene for which accurate light-scattering measure- 
ments of the non-equilibrium temperature fluctuations 
are available [^, • Using information for the thermo- 
dynamic and transport properties [lol 24 1, we have cal- 
culated from Eq. (4) values of pne in the case of liq- 
uid toluene with AT = 10 K at an average temper- 
ature of 25° C. In Table I we compare these values of 
Pne with values for the electromagnetic Casimir pres- 
sure, pem, calculated from Eq. (1), and the critical 
Casimir pressure, Pc = fceTA/L'^ with a Casimir ampli- 
tude A = —0.15 corresponding to periodic boundary con- 
ditions [l^, which are conceptually closest to the stress- 
free boundary conditions adopted in the present paper for 
the non-equilibrium temperature fluctuations. It is seen 
that at L = 1 /um pne is one order of magnitude larger 
than \pc I , and at L = 1 mm is seven orders of magnitude 
larger than |pc|- In terms of a non-equilibrium Casimir 
force per unit area, FneM, with A = 1 mm^, a typi- 
cal value quoted in the literature [1], it means that at 
L = 1 ^m Fne — 1-5 X 1G~* N, and at L = 1 mm 
Fne — 1-5 X 10^^^ N. With modern experimental 
resolutions approaching femto-Newtons [3, [3, US], it 
may even become possible to measure a non-equilibrium 
Casimir effect at a distance L — 1 mm where the more 
traditional Casimir effects are unobservable. Another 
non-equilibrium Casimir effect can be expected in liq- 
uid mixtures where a temperature gradient induces long- 
range concentration fluctuations through the Soret ef- 
fect 0,[2iS| 

In deriving expression (4) for the non-equilibrium 
Casimir pressure, we have not included any possible effect 
of gravity on the non-equilibrium fluctuations. Gravity 
does not affect the intensity of the temperature fluctua- 
tions at small and large wave numbers. However, grav- 
ity suppresses temperature fluctuations at intermediate 
wave numbers in liquid layers with a negative Rayleigh 
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number and enhances temperature fluctuations in liquid 
layers with a positive Rayleigh number 32 . 



This 

means that in a liquid layer heated from above the non- 
equilibrium Casimir pressure could be a bit smaller and 
in liquid layers heated from below larger than that de- 
duced from Eq. (4), as long as the system stays below 
the onset of convection. 
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